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During  the  period  of  the  grant  we  continued  our  studies 
of  one-four  dimensional  classical  and  quantum  lattice  and 
continuous  completely  integrable  systems.  Our  technique  is 
based  on  the  method  of  Backlund  transformations  and  their 
algebraic,  geometric  and  arithmetic  properties.  Methods  of 
Backlund  transformations  were  successfully  used  in  the  study 
of  Pade  approximations.  One  of  our  significant  achievements 
during  the  last  year  was  the  solution  of  the  key  problem  on 
the  almost  perfectness  of  Pad£  approximations  to  solutions 
of  linear  differential  equations. 


I.  Nonlinear  differential  equations  in  dimensions  one. 
two  and  three  with  the  complete  integrability  property. 

The  investigators  studied  three  dimensional  lattice 
and  continuous  models  with  the  complete  integrability  property 
using  Backlund  transformations  (BTs)  algebra  [1],  [2],  [3]# 

[4].  Our  approach,  outlined  in  [1],  is  based  on  topological 
properties  of  symplectic  structures  and  analytic  deformations 
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of  them)  expressed  in  terms  of  commutativity  and  associativity 
axioms  for  S-matrices  associated  with  physical  systems.  This 
topological  approach  is  used  for  classical  and  quantum( opera¬ 
tor  valued)  systems.  Commutativity  and  associativity  axioms 
take  the  form  of  factorization  equations  on  elements  of  S-matrices 
and  are  equivalent  to  algebraic  identities  satisfied  by 
BTs  of  solutions  of  corresponding  nonlinear  differential  equ¬ 
ations.  These  algebraic  identities  satisfied  by  BTs  take  the 
form  of  universal  discrete  completely  integrable  equations 
[1],  [3].  We  present  two  of  these  universal  equations  in 
dimensions  three  and  four  respectively: 


( X_-X^)  _  + 

‘  "  B1  32  '  B3  B' 


(1.1) 


+  ( X^-X  *  n_ 


VBi 


=  Or 


B3  B1  *  B2 


(  X  _i  )  M  )  •  „ 

V  1  4;  V  2  3 '  ~  - 


_  _  * „  +  (X  -X  )  (X  -X  )r_  _  • 

BrB4  B2'B3  VB3  B  2 
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+  (X  -X  )  ( X  -•  )cr_  _  *  _  =  0, 

B1B2  B4B3 


where  is  an  elementary  3T  corresponding  to  the  addition 

Bi 

of  an  apparent  singularity  at  X  =  to  r,  see  [1].  In  a 
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particular  case,  when  completely  integrable  systems  are 

associated  with  second  order  linear  differential  (or  difference) 

spectral  problems,  has  the  form  o_  =  W(y  , .  ..,w )/ 

,  B  B  1  K 

det(X.  .  ,  for  k  eigenfunctions  v.  correspond- 

1  1#3  “  *#*i*#^  1 

ing  to  eigenvalues  X.y  with  BT  data  B  =  (X  ,  ...,X  )y 

X  1  K 

and  a  Wronskian  (continuous  or  discrete)  of  functions 

T^,...,^.  Discrete  equations  (1,1),  (1.2)  take  the  form  of 

the  law  of  addition,  when  BTs  are  represented  as  translation 

operators  in  the  auxiliary  space  of  infinitely  many  variables 

x^:  n  =  1,2,...  (being  Newton's  symmetric  functions),  with 

BT  action  cr  def)  exp  {  E*  X n  ®  or 

i  n 

o_  (X)  =  c(X  +  X.).  The  equations  (1.1),  (1.2)  in  this 

B.  1 

i 

form  are  equivalent  to  the  laws  of  addition  on  curves  (of 
genus  g^T  "),  and  are  crucial  in  the  complete  solution  of  the 
Schottky  problem  of  the  determination  of  Jacobian  varieties 
among  all  Abelian  varieties. 

Discrete  universal  equations  (1.1)  generate  in  various 
limits  and  reductions  the  Kadomtzev-Petviashvili  (KP)  equation 
and  various  other  well  known  two  dimensional  completely  inte¬ 
grable  systems,  see  [5].  We  have  proved  in  [1],  [4]  that 

all  multicomponent  (operator)  two  dimensional  completely 
integrable  systems  of  isospectral  deformation  nature  are 
algebraic  reductions  of  the  universal  equation  (1.1).  Moreover, 
particular  algebraic  reductions  of  (1.1)  that  determined  two- 
dimensional  multicomponent  completely  integrable  systems  were 
completely  described  in  [4]  using  infinite  dimensional  algebras 
arising  from  particular  algebras  of  pseudodifferential  opera¬ 
tors.  Recently  we  were  able  to  prove  that  any  multicomponent 
fnatrix,  operator)  three  dimensional  systems  with  the  complete 
integrability  property  also  arise  from  systems  (1.1)  with 


a  particular  representation  of  the  algebra  of  BTs.  These 
three  dimensional  matrix  systems  include  any  system  of  non¬ 
linear  p.d.e.  that  arises  as  a  commutativity  condition 


[L  ,  L2]  =  0,  where 


L^,  L2  are  linear  differential  operators 


in 


&x 


bz 


efficients  (and  are  of  the  first  order  in 


with  matrix  (operator)  co- 
l  b 


)  .  Thess 


3z 

three-dimensional  matrix  systems  include  the  general  monopole 
equations  for  SU(n)  or  S0(n)  cases  and,  say,  the  matrix 
generalizations  of  the  KP  equation  derived  by  the  authors  in 

[6]. 


The  investigators,  together  with  M.  Tabor,  succeeded 
in  proving  the  "Painleve  property"  for  two  dimensional 
multicomponent  isospectral  deformation  equations  including 
various  versions  of  the  matrix  nonlinear  Schrodinger, 
mdV  and  KdV  equations  [1] .  This  led  to  the  establish¬ 
ment  of  new  completely  integrable  many  particle  systems 
arising  in  the  pole  expansion  of  meromorphic  solutions  (as 
a  particular  case  of  Painleve* s  expansion).  These  results 
were  based  on  the  above  mentioned  reductions  to  the  univer¬ 
sal  equation  [4],  [7] .  The  investigation  of  the  Painleve 

property  is  extended  to  matrix  three  dimensional  systems. 
Simultaneously,  we  arrived  at  a  solution  of  a  new  algebraic 
problem  of  the  determination  of  all  exponents  of  various 
branches  of  the  Painleve  expansions  of  solutions  of  various 
flows  commiting  with  a  given  completely  integrable  one 
(already  a  nontrivial  problem  for  the  n-th  KdV  flow) . 

Painleve  equations  without  movable  critical  points 
were  found  by  the  investigators  to  be  of  utmost  importance 


II.  Pade  tyoe  approximation  methods  and  their  appli¬ 
cations  in  mathematics  and  mathematical  physics. 

We  are  continuing  to  be  active  in  theoretical  and 

computer  studies  of  Pade  approximations  to  solutions  of 

differential  equations.  Here,  under  the  generic  name  of  Pade- 

type  approximations,  we  understand  the  approximation  of  functions 

having  convergent  or  asymptotic  power  series  expansions, 

locally  by  rational  functions.  If  we  consider  approximations 

in  only  nonarchimedian  metrics  (  |t  I  =  pi™  iff  f(x)  = 

.Co  1 

E  an(x  -  xQ)n  and  m  =  min  tn  :  an  ^  01)  .  Then  we  obtain 
classical  (multipoint)  Pade  approximations.  Our  most  spec¬ 
tacular  results  during  the  last  year  concern  Pade^-type 
approximations  to  solutions  of  linear  differential  equations. 
Pade-type  approximations  in  this  case  are  directly  connected 
with  the  algebraic  formalism  of  the  inverse  scattering 
method  and  Backlund  transformations  (2],  [3],  [9],  [10].  Pade 

approximations  to  solutions  of  linear  differential  equations 
are  very  attractive  in  view  of  the  possibility  of  establishing 
rigorously  the  convergence  of  this  method  for  multivalued 
solutions  of  differential  equations  following  the  monodromy 
properties.  The  efficiency  of  Pade  and  Pade-type  approximations 
is  also  demonstrated  by  the  fact  that  these  methods  furnish 
exceptionally  good  rational  approximations  to  numbers  that  are 
values  of  solutions  of  differential  equations.  Sometimes  we 
obtain  explicit  expressions  of  the  continued  fraction  expan¬ 
sions  of  these  numbers.  For  example,  we  studied  new,  Pade- 
type  approximations  to  such  functions  as  ln(l  -  1/z)  at 


arctg 


In  this  way  v:2  obtain 


explicit  expressions  for  the  best  rational  approximations  to 

TT 

transcendental  numbers  such  as  ln2  and  ^  .  See  the 

Appendix  for  explicit  three-term  recurrences  relating  rational 

TT 

approximations  to  .  These  rational  approximations  give 

outstanding  bounds  for  measures  of  diophantine  approximations 
to  classical  transcendental  numbers  [8],  [11]#  [12]. 


Our  main  result  [14],  [15]  on  the  theory  of  Pade* 

approximations  itself  is  the  proof  of  the  "almost  almost  perfectness' 

of  Pade  approximations  to  solutions  of  arbitrary  (algebraic) 

differential  equations.  Our  result  [14],  [15]  gives  a  complete 

solution  to  Kolchin's  problem  [13]  that  for  any  solution 

£(x)  of  an  arbitrary (algebraic)  differential  equation  over 

C(x)  regular  at  x  =  any  e  >0  and  any  arbitrary  rational 

function  P(x)/Q(x),  we  have  ord  __  ff(x)  -  P^X]/q{x)J 

x  -  *  »  ' 

^  (2+5)  max  (deg(P),  deg(Q)}  +  CQ(5,f).  We  want  to  note 
that  C.F.  Osgood  recently  announced  that  he  too  has  an 
effective  solution  to  Kolchin's  problems  in  the  case  of  solutions 
of  linear  differential  equations. 

The  general  theorem  [14],  [15]  of  the  investigators 

shows  that  Pade  approximations  are  "almost  almost  perfect." 

To  invoke  the  notion  of  perfectness  we  remind  of  the  struc¬ 
ture  of  the  Pad£  table  for  f(x).  On  the  (n,m)th  place  of 

the  Pade  table  for  f(x)  at  x  =  x  ,  we  put  a  rational  function 

o 

Pn(x)/ >  ,  n’  such  that  ord  x=x  (f(x)  -  P  fx)/Q  (x>)  >  n+m. 
n  Qm(x)  0  n  m 

Then  the  Pad£  table  is  divided  into  squares  consisting 

of  identical  rational  functions.  The  standard 


*  ,v 


( 


definition  of  perfectness (normality)  means  that  each  square 
is  of  size  one.  Unfortunately,  this  natural  definition  is 
not  universally  valid  (e.g.  it  is  enough  to  multiply  any 
f(x)  by  an  appropriate  rational  function  to  create  squares  of 
an  arbitrarily  large  size) .  A  more  reasonable  condition  is 
almost  perfectness  (Mahler,  1935  [16],  [8])  meaning  that 
squares  in  the  Pad£  table  have  uniformly  bounded  sizes.  This 
is  equivalent  to  e  =  0  in  the  solution  of  the  Kolchin  problem 
above,  while  the  investigators'  solution  [14],  [15]  of  the 

Kolchin  problem  is  only  an  "almost  almost  perfectness" 
statement.  The  investigators  now  have  the  proof  of  almost 
perfectness  for  a  large  class  of  linear  differential  equations, 
including  differential  equations  of  hypergeometric  type,  and 
for  arbitrary  algebraic  functions.  The  only  analytic  results 
that  existed  before  were  those  of  Arms  and  Edrei,  see  [17], 
for  particular  trigonomethric  functions  such  as  cos  x,  —  x 

using  the  positivity  argument. 


Our  results  are  generalized  to  Hermite-Pade  approxima¬ 
tions  to  sequences  of  functions  [14],  [15].  These  results 
and  the  methods  of  algebraic  geometry  and  differential  algebra 
that  we  use,  were  already  applied  by  the  investigators  [15], 

[18]  to  study  the  diophantine  approximations  of  values  of 
solutions  of  differential  equations.  In  particular,  as  a 
solution  of  Lang's  problem  we  proved  the  best  possible  "2  +  e" 
bound  in  the  measure  of  irrationality  of  any  value  9  of 
E-function  f(z)  at  rational  z  =  r.  Here  the  "2  +  «" 

bound  means  that  for  any  e-  >  0,  [  ?  -  P/q  J  >  |q 

for  rational  integers  p,q  with  |q|.*q  (3,  e)  An  E-function 

a  a 

f(z)  is  a  function  f(z)  =  ”-  -r-  xn.  with  a  6  0, 

n=o  n I  n 


satisfying  a  linear  differential  equation  with  rational  func¬ 
tion  coefficients,  and  9  =  f(r),  r  ^  0. 


Prom  the  point  of  view  of  rational  approximations, 
our  results  for  e  =  0  imply  the  boundedness  of  (degrees  of) 
elements  in  the  continued  fraction  expansion  of  solutions  of 
linear  differential  equations.  These  results  and  the  solution 
of  the  Kolchin  problem  are  the  first  step  in  the  explicit 
determination  of  rates  and  domains  of  convergence  of  rational 
approximations  to  solutions  of  arbitrary  differential  equa¬ 
tions  . 


Our  work  was  facilitated  by  use  of  symbolic  algebra 
manipulation  systems  SCRATHPAD  at  IBM  and  one  of  the  versions 
SM?  (Caltech) . 
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Appendix 


We  present  the  explicit  expressions  for  the  "good" 

4rr 

rational  approximations  to  ,  and  the  three-term 

recurrence  that  these  approximations  satisfy.  The  explicit 
formulas  are  specializations  of  Pad^-type  approximations  to 


8  1 
rt  arot9  7r-- 

tions  X  /Y  to 
n  n 


at  z  =  3.  In  the  rational  approxima¬ 


tions  X  /Y  to  ,  both  sequences  X  and  Y  are 

n  n  313  ^  n  n 

solutions  of  the  following  three-term  recurrence: 


(2.1) 


Vn>  Zn+2  +  V">  Zn+1  +  Vn)  Zn  '  ° 


for  Zn  =  X^  or  =  Y^:  n  =  0,  1,  2,...  .  The  coefficients 
AQ(n),  A^(n),  A2(n)  are  Po3-ynomials  of  degree  9  in  n: 


A2(n)  *  -8 *  (4n  +  7)  •  (4n  +  5)*(4n  +  3)  •  (4n  +  l)'(2n  +  3) 
.  (n  +  2) * ( 27 279n3  +  52l64n2  +  31511n  +  6946); 


AL(n)  =  3  •  (4n  +  3) *  (4n  +  1) •( 154846 24231n  + 


& 


1225 1806648 2n6  +  401859 218160n5  + 


706125904254n4  +  715 28 2318379n3  + 
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we  obtain 


Using  the  approximation  X  /Y  to  - — _ 

n  n  3lT 

the  following  bound  for  the  measure  of  irrationality  of 

:  I71-/  / 3  -  p/q  [  >  |q  |  4-81‘'  !  for  all  rational 

integers  p,q?  |q  I  >  qQ. 


